DISPERSION FORCES AND DUALITY 



STEFAN YOSHI BUHMANN, STEFAN SCHEEL 

Quantum Optics and Laser Science, Blackett Laboratory, Imperial College London, Prince 
Consort Road, London SW7 2AZ, United Kingdom 

s.buhmann&imperial. ac.uk 

HASSAN SAFARI, DIRK-GUNNAR WELSCH 

Theoretisch-Physikalisches Institut, Friedrich- Schiller- Universitdt Jena, Max-Wien-Platz 1, 

07743 Jena, Germany 

Received November 26, 2008 

We formulate a symmetry principle on the basis of the duality of electric and magnetic 
fields and apply it to dispersion forces. Within the context of macroscopic quantum 
electrodynamics, we rigorously establish duality invariance for the free electromagnetic 
field in the presence of causal magnetoelectrics. Dispersion forces are given in terms 
of the Green tensor for the electromagnetic field and the atomic response functions. 
After discussing the behavior of the Green tensor under a duality transformation, we are 
able to show that Casimir forces on bodies in free space as well as local-field corrected 
Casimir-Polder and van der Waals forces are duality invariant. 
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1. Introduction 

Dispersion forces suclt as tlte Casimir forc^ on a body, the Casimir-Polder (CP) 
forcJ2l between an atom and a body and the van der Waals (vdW) forcJ^I between 
two atoms were originally conceived as effective electromagnetic forces between elec- 
trically neutral, but polarizable ground-state objects; they are typically attractive.^ 
Somewhat later, the investigations were extended to bodies and atoms with an addi- 
tional magnetic response^ revealing that polarizable and magnetizable objects may 
repel each other. Forces between magnetoelectric systems have recently been subject 
to a renewed interest^ due to the availability of metamaterials with a controllable 
permittivity and permeability.^^ 

As we will demonstrate in this article, investigations of this kind can be con- 
siderably simplified by exploiting the well-known duality of electric and magnetic 
fields.'^ To that end, we first study the effect of duality transformations in the con- 
text of macroscopic quantum electrodynamics*^- (Sec. [21) and then use our results to 
prove that duality invariance is a valid symmetry of dispersion forces under very 
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general conditions (Sec. [3]). 



2. Macroscopic quantum electrodynamics and duality 

Duality is one of the inherent symmetries of the Maxwell equations. To see this, 
we group the electric and magnetic fields E{r,t), B{r,t) and excitations D{r,t), 
H{r,t) into dual pairs {E, ZqH)'^ and {ZqD, B)"^ , where the vacuum impedance 
Zq = v^/io /eo has been introduced for dimensional reasons. In this dual-pair nota- 
tion, the Maxwell equations in the absence of free charges and currents assume the 
compact form 



Grouping the polarization P{r,t) and the magnetization M{r,t) according to 
{ZqP, fioM)"^ , the relation between the fields and excitations reads 



ZoD\ ^If E \ f ZoP\ 
B ) c\ZoHj UoM/ 



(3) 



It is now immediately obvious that Eqs. (H])-® are invariant with respect to a 
duality transformation 

^Vir,9)(^), Vir,e)^( ^'^"^^^^^^J G x S0(2), (4) 

with I?(r, 9) being the most general real matrix that commutes with the symplectic 
matrix in Eq. This transformation may be viewed as a rotation in the space 
of dual pairs (0 < 6* < 2tt) together with a rescaling of all fields (r > 0). In classi- 
cal physics, the real-valued electromagnetic fields are often combined into complex 
Riemann-Silberstein vectors x+iy, in which case duality invariance manifests itself 
as a U(l) symmetry. 

Let us next address the compatibility of the duality transformation with the 
constitutive relations, which may conveniently be formulated in terms of the Fourier 
components x of the fields, x{r, t) — duo x[r , uj, t)+h..c. For linear, local, isotropic, 
dispersing and absorbing media, the constitutive relations may be given as 



Zob\ ^lfeO\f E \ flO\f ZoP^ 
B J c\Ofi)\ZoHj ^Om/UoMn 



(5) 



where e = e(r,w) and ^ — ii{r,uj) denote the relative electric permittivity and 
magnetic permeability; and and Mj^ are the noise polarization and magnetiza- 
tion which necessarily arise in the presence of absorbing media. Invariance of the 
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constitutive relations under the duality transformation requires that 

^* ^\ = V(r f)^f^ ^^-n-^fr '^cos^e + ^sin^e (/x - e) sin 61 coi 

) ^'\OfiJ UAi-£)sinecos0 esin^^ + ^cos' 



\ ( rcos6 fir sin 9 \ / P^j^ 

M^/c) " l^-(l//^*)rsin0 (fi/ fi*)r cos 9 ) \M^/c 



(6) 
(7) 



Condition ^ can be fulfilled in two ways: ft holds if the relative impedance of the 
media is equal to unity, Z = \J~[sJs= 1. In this case, which includes both free space 
and a perfect lens medium (£ = ^ = — 1}^, the duality rotations form a continuous 
SO (2) symmetry of the electromagnetic field and one has e* = [i* — e as well as 



Pj^ \ ( rcos9 ersin6'\ / P^ 



M_^/cJ \-{lle)rsin9 rcos9 J\M_^/c 



(8) 



For media with a nontrivial impedance, Eq. ([6]) holds for 9 = 1111/2 with ?igZ only. 
The presence of such media hence reduces the duality invariance from the full SO (2) 
group to a discrete Z4 symmetry with the four distinct members 

V^=I, = ( J) , = -I, = -Pi, (9) 
(X: unit matrix) where Eqs. ([5]) and ([7]) imply the transformations 



£ \ / COS^ 9 SIT? U \ I £ 



/i y \ sm^ 9 cos^ tl J \ii 

Pf^ \* / rcos9 ^rsin9\ / 
Mj^/c) ^ \-{l/£)rsm9 rcos9 )\M_^/c 



(10) 
(11) 



Duality is thus an exact symmetry of the Maxwell equations in the absence of 
free charges and currents. It must also be manifest in the underlying Hamiltoniarff^ 
Hp = J2\=e m I 'l^^ Iq° '^^ ^ fx ' /a ; whcrc the bosonic dynamical vari- 
ables fe(r,uj), fm{r,uj) are associated with electric and magnetic medium-field 
excitations. Noting that the dynamical variables are related to the noise fields via 



UoMi,; V ^ ^ vwi/iMiyv/n. 



(12) 



and recalling Eqs. ([8]) and (fTT|) . they are seen to transform as 

' ^ (13) 



fe\ ^ f rcos9 -i{n/\n\)rsin9\/ fe 
fra) V-i(|£|/e)rsin6l rcos6' )\frn 

in both the continuous {e = ^) and discrete (6' = n7r/2) cases. It follows that a duality 
transformation leads to a rescaling of the Hamiltonian Hp = r^Hp, such that the 
equations of motion remain invariant. 
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It is important to note that electromagnetic forces are not duality-invariant in 
general, even when acting on electrically neutral systems. For instance, the Lorentz 
force on a neutral magnetoelectric body of volume V can be written aJ^l 



/ dA- {eoE{r)E{r) + —B{r)B{r) - \ 
Jdv I A^o ^ 



Mo 

-eo^y" dV£;(r)xB(r) (14) 

(I: unit tensor); it is obviously not duality-invariant. Duality invariancc would be 
realized for a stationary field acting on a body at rest (such that the total time 
derivative vanishes), provided that D eqE and HfuB/fiQ on the body surface. 
While this can never be true on an operator level due to the unavoidable presence of 
the noise fields and AiN, one would expect the Casimir force (F) on a stationary 
body in free space to be invariant. The situation is very similar for the Lorentz force 
on a neut ral at om (position r^, polarization Pa, magnetization Ma) which can be 
written aJSM 

F = Va J dV PAir)-E{r) + MAir)-B{r) + PAir)xi'A-B{r) 

+ d\PA{r)xB{r) (15) 

when neglecting diamagnetic interactions. Noting that Pa and Ma transform under 
duality like P and M, one sees that duality invariance can only hold for an atom 
at rest (so that the velocity-dependent terms do not contribute) prepared in an 
incoherent superposition of energy eigenstates and subject to a stationary field 
(so that the time-derivative does not contributill|), provided that D « SqE and 
H « B / fia hold within the volume occupied by the atom. Again, this is never 
true on an operator level, but one may expect CP and vdW forces on atoms to be 
duality invariant. We will confirm the conjectured duality invariance of dispersion 
forces in Sec. |3l In the following, we concentrate our attention on the particular 
transformation r = l, 9^t:/2 which is a generator of the discrete duality group Z4. 



3. Duality invariance of dispersion forces 

As a preparation for studying the duality invariance of dispersion forces, let us first 
establish a few general expressions for these forces in terms of the relevant response 
functions. By taking the ground-state expectation value of the Lorentz force (fT4|) . 
one easily finds that the Casimir force on a stationary homogeneous magnetoelectric 
body is given bjff^l 



//•oo / 2 



Tr<^ I X 



Vx Vx 



(16) 
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where G^^'' is the scattering part of the Green tensor G of the electromagnetic field, 
1 



Vx 



UJ 



Vx J e{r,Lu) 



G(r, r', Lj) = 6{r — r'). 



Alternatively, the Casimir force may be given as a surface integraff^l F^F^ 

TT 



(17) 







dA- 



dV 



;(!)( 



r,r,iO-ilTrGi\)(r,r,iO (A = e,m) (18) 



with Gee(r, r', w) = (iw/c)G(r, r', w)(iw/c), Gmm(r, r', w) = V x G(r, r', w) x V. 

The CP force on a single atom can be derived from the Casimir force in its 
volume-integral form p6p by considering the force on a dilute gas of atoms [num- 
ber density r/(r), polarizability a(a;), magnetizability /3(w)] occupying an otherwise 
empty volume V . Within leading order in 77, the contributions of the atoms to 
the permittivity and the inverse permeability {n — 1^^^) is given by the linearized 
Clausius-Mosotti relationJI^I Agf r, u)) — ri{r)a{u)) / £n and ls.n{r ,uo) — —ri{r)P{uj)^,Q. 
Using a linear Born expansion^'^^l^ one finds that the resulting change of the Green 
tensor reads 

AG(r,r',u;) = 



3 ?7(s) iuj 



£0 



■ a{uj)G{r, s,uj)-G{s,r' ,uj) 



G{r,s,uj)xVs ■ VsxG(s,r',w) 



(19) 



Combining this with Eg. ([TC)) one can show that the Casimir force on the atomic 
cloud can be written aa^i^ = ^ /y 7](r) VC/ (r) where U {rA) — Ue{rA)+Ue{rA) 
with 



27r£o 



/ d^aA(iOTrGi\^(r^,r^,iO (A = e,m) 
Jo 



(20) 



{ae = ct, ctm — P/c^) is the CP potential sought. The vdW potential between two 
atoms A and B can be obtained in an analogous way by introducing a second 
dilute cloud of atoms and applying Eq. (fT9|) to the CP potential. This results in 
U{rA) = d'^r 77(r)[/(r^, r), where the required two-atom vdW potential reads 

U{rA-,rB) = Uee{rA,rB) + Uern {v A,r b) + Ur,ie{r A-, T b) + Urnm{r A,r b) , 



U\x'{rA,rB) 



'0 JO 



deaA(iO«f'(iOTr 



Gaa' (t-a, t-b, iC) •GA'A(r"B, n4, iS.) 

(A,A' = e,m) (21) 



with Gem{r, r', w) = {iuj/c)G{r, r', uj) x V, G™e(r, r', lj) = V x G(r, r', w)(icj/c). 

Dispersion forces and potentials can thus be given in terms of the response 
functions of the electromagnetic field and the atoms, so their behavior under a 
duality transformation can be determined from that of G, a and /3. By virtue of the 
linearized Clausius-Mosotti relations, the known transformation properties e* = 
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fjL* = e imply that a* = (3* — a/c^. As shown m |Appcndix A[ the Green tensor 
transforms according to 

G:,(r, r', u) = ^i-\r, c^)G„„n(r, r', LJ)^l-\r' , uj) + ^i-\r, uj)Sir-r'), (22) 



r', uj) = e(r, uj)Gee{r, r' , uj) e(r', uj) - £(r, uj)S{r-r'), 
Ge„(r,r',w) = -^"i(r,w)G„e(r,r',w)£(r',w), 
GrneC'^;'''.^) = -e{r,Lj)Gem{r,r',uj)fi^'^{r',u}). 



(23) 
(24) 
(25) 



These laws immediately show that the discrete global duality transformation 
s ^ fi, a ^ fi/c? leaves ground-state dispersion forces on stationary objects (and 
associated potentials) invariant, where the individual electric and magnetic compo- 
nents ([T51) , (HOI) and (PT|) transform into one another according to ^ Fm , Ue Um 
and Uee^Umm, Uem^Ume, respectively. Due to the factors e and /i^^ appearing 



in Eqs. l|22p - (P5l) . this only holds for forces on atoms and bodies which are situated 
in free space. 

In order to extend duality invariance to atoms which are embedded in a medium, 
local-field corrections need to be taken into account. Using the real-cavity model, 
one can show that local field effects give rise to correction factorJl^ 



Ce{uj) 



1 2 



2e{Lj) + 1 

so the potentials ^U\i and (I^T]) generalize to 

h 



2/i(cj) + 1 



UxivA) 

Uw(rA,rB) 



27reo 

h 



27re2 



/ decA(iOaA(ie)TrGW(r^,r^,iO 
Jo 

dec^(iOcf'(iO«A(iO«f'(ie) 



(26) 



(A = e,TO) (27) 



X Tr 



"'0 



'AA 



(rA,rB,iO-GA'A(rB,r-yi,iO 7 (A,A' = e,m). (28) 



When applying a duality transformation, the factors e and /k"^ arising from the 
transformation of the Green tensor ((22l) - ((25)) combine with those contained in the 
local-field correction factors ((26)l in such a way that the corrected potentials ((27|) and 
(|28p transform into one another according to Ue Um, Uee ^ Umm and Uem ^ Ume- 
When including local-field corrections, the total ground-state dispersion potentials 
are hence also duality invariant for embedded atoms. 



4. Summary 

We have studied the behavior of electromagnetic fields, response functions and dis- 
persion forces under duality transformations. In the presence of media with a non- 
trivial impedance, the SO (2) symmetry of duality rotations reduces to a discrete Z4 
symmetry. The duality transformations induced by the generator of this group are 
displayed in Tab. [TJ Note that the transformation needs to be applied four times 
to the fields (first block) in order to return to the original state. On the contrary. 
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Table 1. Effect of the duality transformation with r = l,S = 7r/2on fields, response 
functions and dispersion forces. 



Dual partners 


Transformation 




E, H: 


E = ZqH, 


H* = -E/Zo 


£>, B: 


D* = BZo, 


B* = -ZoD 


P, M: 


P* = M/c, 


M* = -cP 


Pa, Ma- 


Pi = Ma/c, 


MX = -cPa 




_ PfJ = (iMt^/c^ 


M* = -cP^/e 


fe , fm '- 


ft = -iWI^D/m, 


/* = -i{\e\ls)h 


e, fi: 


e* = fi, 


11* = e 


a, /3: 


a* = c^fi. 


13* = a/c^ 


Gee ) Gmm ■ 


gl = (i/m)G„™(i/m) + (i/m)5, 


Gmm = eGeeE - sS 


Gem; Gme- 


GJm = -(l//i)Gmee, 


Gme = -eGem(l/At) 




F* = Fm, 


Ft^ = Fe 


Ue, Um- 


U: = Um, 


UJn = Ue 


Uee 1 Umm '■ 




^mm — 






U-me — Uem 



the transformation is self-inverse when acting on the response functions of bodies, 
atoms and the electromagnetic field (second block). We have shown that the electric 
and magnetic components of dispersion forces on neutral bodies and atoms, which 
are at rest and situated in free space, depend on these response functions in such a 
way that they inherit a similar transformation behavior (third block). As a conse- 
quence, the total dispersion forces are duality invariant. As demonstrated, duality 
invariance can be extended to atoms embedded in media provided that local-field 
corrections are taken into account. 

The duality invariance established in our work provides an important consistency 
check for investigations of dispersion forces. It can further serve as calculational tool: 
Once the force on an object in a particular magnetoelectric environment is known, 
an expression for the force in the dual arrangement can be generated by simply 
making the replacements a <-> fijf? , £ <-> /i. 
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Appendix A. Duality transformation of the Green tensor 

To derive the transformed Green tensor G*, which is a solution to the Eq. (fT7|) with 
e* = n and ji* =e instead of e and /i, we first note that the Maxwell equations ([1]) 
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and (HJ together with the constitutive relations (O are uniquely solved by 



-1 /"dV'Gee(r, r', u) -P^ir', uj)^ZoJ dV'G„„(r, r', lo) -M^ir', lu), 
'° (A.l) 

-Zo I d^r'G,ne{r,r',oj)-Pj^{r',Lj) - fio / d\'Gm^{r,r' ,u;)-Mj^{r' ,uj), 

(A.2) 

e{r,uj) 



dV 



dV 



e{r, uj)Gee{r, r' , us) — 6{r — r') 
fd\'G^e{r,r',Lu)-P^{r',uj) 



Gm,n{r,r',uj) 



+ S{r-r') 



(A.3) 



(A.4) 



Applying the duality transformation to Eqs. (jA.lD and (jA.2p with the aid of the 
transformation laws established in Sec. [Hand using Eqs. (|A.3p and (jA.4[) . the un- 
known quantities G*xx' on the rhs of the transformed equations can be related to 
the untransformed ones G\y appearing on the Ihs, and one obtains Eqs. 
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